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Abstract 

Hartle’s model for slowly rotating stars has been extensively used to compute 
equilibrium configurations of slowly rotating stars to second order in perturbation 
theory in General Relativity, given a barotropic equation of state (EOS). A recent 
study based on the modern theory of perturbed matchings show that the model must 
be amended to accommodate EOS’s in which the energy density does not vanish at 
the surface of the non rotating star. In particular, the expression for the change in 
mass given in the original model, i.e. a contribution to the mass that arises when 
the perturbations are chosen so that the pressure of the rotating and non rotating 
configurations agree, must be modified with an additional term. In this paper, the 
amended change in mass is calculated for the case of constant density stars. 


1 Introduction 

Hartle’s model [6] describes the axially symmetric equilibrium configuration of a rotating 
isolated compact body and its vacuum exterior in perturbation theory in GR. The interior 
of the body is a perfect fluid satisfying a barotropic equation of state, it does not have 
convective motions and rotates rigidly. The exterior is an asymptotically flat, stationary 
and axisymmetric vacuum and both the interior and the exterior share equatorial symme¬ 
try. The interior and exterior are matched across a timelike hypersurface that represents 
the boundary of the star. 

Hartle’s model is built taking slow rotation perturbations around a static and spher¬ 
ically symmetric background. The perturbed field equations for the interior region are 
presented with conditions at the origin so that, once a barotropic EOS is specified, reg¬ 
ular solutions for the fluid are found. On the other hand, the vacuum equations can be 
analytically solved and the solutions compatible with asymptotic flatness are given. The 
solutions for both fluid and vacuum are completely determined assuming that the metric 
perturbations are continuous in the coordinates used at the boundary of the star, where 
the fluid and vacuum spacetimes are matched. 

However, the theory that provides the matching conditions in a perturbed sense to 
first and second order, formulated independently of the coordinates used to describe the 
background configurations and also of the gauges (both spacetime and hypersurface) 
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used to describe the perturbed configurations, is very recent [8]. Hartle’s model has 
been revisited in a fully consistent manner from the perspective of the modern theory of 
perturbed matchings in ms- In particular, in that work it has been found that a function 
in the second order perturbations, m 0 (r), in the coordinates and gauges used in the original 
work |[6j is not continuous, in contradiction with the implicit assumption there. It turns 
out that the discontinuity of the function mo(r) affects a particular outcome of the model, 
the change in mass 5M, which is the contribution to the mass of the rotating configuration 
that arises when the central densities (pressures in this paper) of the non-rotating and 
rotating configurations are forced to agree. The jump on a radial m 0 function already 
appears implicitly in the revision of Hartle’s model in [1], where the matching hypersurface 
is prescribed as the set of points where the perturbed pressure vanishes in a certain gauge. 

The correction to 5M found in [10] is proportional to the energy density evaluated at 
the (non rotating) surface of the star. Hence, it is negligible when the usual equations 
of state for neutron stars are considered, since the pressure and the energy density typi¬ 
cally decrease together and vanish simultaneously at the surface of the star. This is the 
behaviour shown by, e.g. polytropic EOS’s. Nevertheless, the correction may be relevant 
in other EOS’s for which the energy density takes a finite value at the boundary, which is 
precisely the case of linear equations of state, for instance those used to describe strange 
quark stars [3], or constant density (homogeneous) stars [3]. 

In this work the particular case of homogeneous stars is studied. This equation of state 
may not be realistic in physical terms, but its study is interesting for several reasons. First 
of all the spherically symmetric and static configuration can be solved analytically and 
hence, the perturbed field equations become simpler. 

Secondly, it is interesting to know how noticeable the correction for the change in mass 
might be in numerical terms. There are factors in the change in mass that are easy to 
estimate and, in fact, can be taken as inputs for the model, as the mass and the size of 
the static and spherically symmetric star. In contrast, the perturbation to the pressure is 
not easily estimated and the model must be solved to second order. Another important 
factor is the value of the energy density at the surface of the star in the nonrotating 
configuration and this is precisely the reason for having chosen this particular EOS. It is 
probably one of the most favourable cases for the correction. Thus, it is reasonable to 
think that the constant energy EOS may constitute a numerical bound for the amended 
change in mass, since for any other of the usual EOS the value of the energy density at 
the surface will not be as important as in the present case. 

Homogeneous stars drew the attention of Chandrasekhar and Miller [3], back in 1974. 
In that work, they used Hartle’s formalism to solve the homogeneous star in the slow 
rotating approximation. These have also been studied under the CMMR formalism [2], a 
treatment based on the post-Minkowskian and small deformation approximations in [5], 
as a particular case of the linear EOS. 

In this paper Chandrasekhar and Miller’s work is retaken, following Hartle’s setting 
as described in ra and taking into account the correct value of 8M. The results are 
presented in an analogous way to [3J to ease the comparison of tables and figures. 
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2 Hartle’s model in brief 

Hartle’s model is based upon the following one parameter family of metrics 0 

g £ = — e v ^ (l + 2 e 1 2 h{r, 9)) dt 2 + e A( ^ ^1 + 2e 2 —— m(r, 9)^j dr 2 

+r 2 (l + 2 e 2 k(r, 9 )) [d9 2 + sin 2 9(dip — eu;(r, 9)dt) 2 ~\ + 0(e 3 ), 

defined on the static and spherically symmetric background spacetime (Vo ,g := g e =o) so 
that 

g = -e^dt 2 + e A(r W 2 + r 2 (d9 2 + sin 2 9dy 2 ). (1) 

The first and second order metric perturbation tensors, K\ = d e g e | 0 and K 2 = 
d 2 g e \ £=0 take the form 

Ki = — 2r 2 u{r) sin 2 9dtd(p : (2) 

. , 4 p 2 A(r) 

K 2 = (—4 e u ^h(r, 9) + 2r 2 sin 2 9u 2 (r , 6*)) dt 2 H---m(r, 6*)dr 2 

+4r 2 fc(r, 9)(d9 2 + sin 2 9dip 2 ), (3) 

where the functions in ft 2 admit the decompositions h(r,9) = h 0 (r) + h 2 {r)P 2 {cos 9), 
m(r,9) = m 0 (r) + m 2 (r)P 2 (cos 9) and k(r,9) = k 2 (r)P 2 (cos 9). Note that the choice 
ko{r) = 0 fixes the second order spacetime gauge at both regions (fluid and vacuum, 
see [ID]). The second order perturbation tensor as written in (|3]) with fco(r) = 0 is said 
to be in the fc-gauge in [TO]. The perturbed Einstein’s field equations are imposed as 
follows. Consider the e family of equations G(g £ ) a p = 87rT £a p and take the first order 
field equations as the first e-derivative at e — 0, i.e. d £ G(g £ ) a p\ £= Q = d £ T £a8 \ £=Q and the 
same for successive orders. For the interior region consider a family of perfect fluids so 
that T e = (E e + P £ )u £ ®u £ + P £ g e . The unit fluid flow under the assumption of circularity 
and rigid rotation explicitly reads 

u £ = e v ^ 2 dt + + —e 3v ^ 2 {VL 2 g w + + K 2 tt/‘2)dt + 0(e 3 ), (4) 

for some constant Q. The energy density and pressure admit the expansions 

2 2 

E e (r,6) = E{r ) + jE^ 2 \r,9), P £ (r : 9) = P(r) + yP (2) M), (5) 

where the first order terms in e vanish due to the structure of the perturbed field equations 
to first order [EjJ (see also [U [TO]). The field equations lead to E^ 2 \r,9) = E^\r) + 
E 2 2 \r)P 2 (cos9 ) and P (2 ^(r, 9) = Pq 2 \t) + p!^\r)P 2 {cos9) [6] (see also [T ID]). 

Consider two copies of the spacetime introduced so far, (V^, g + ) with the correspond¬ 
ing family gf to describe the perfect fluid region, and (Vj - ,# - ) with g~ for the vacuum 
region. The coordinates are labelled with subscripts in order to avoid confusion. 

1 The function m{r,8) in (m) is the same function m(r, 9) that is used in [Bj. In contrast, the m(r,9 ) 

in da corresponds to the re*( r hn(r,9) here. 
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The background configuration for the interior region satisfies the equations of general 
relativistic hydrostatics addressed in [6] 


dM + (r +) 
dr .|_ 

dP(r + ) 

dr + 


4irr + 2 E(r + ), 

(E(r+) + P(r+))(M+(r+) + 4 tt r + 3 P) 
r + (r + — 2M+(r + )) 


( 6 ) 

(7) 


where the functions M + (r + ) and P(r + ) are related to the metric functions A + (r + ) and 
z/ + (r + ) by 


e -A + (r + ) = j_2M>+) 

' r + 

dv + (r + ) 2 dP(r + ) 

dr .|_ P(r + ) + P( r +) dr + 


( 8 ) 

(9) 


This set of equations determines the interior configuration once the values of the energy 
density E and pressure P at the origin are given. The function v + is determined up to 
an additive constant. 

For the exterior vacuum E = P = 0, thus M - (r_) — M is a constant and the metric 
functions are given by 


e "-(r-) _ e -Mr-) = 1 _ 2M 


( 10 ) 


These two background spacetimes (V^,# 4 ") and are matched across diffeo- 

morphic boundaries, So := E^. The most general matching preserving the symmetries 
JllJ (spherical and static) can be cast in parametric form in terms of three coordinates 
{t, (ft, $} on E 0 as 

E o = (4 = T ,<P+ = r + = R +, °+ = i s o = {*- = T , ( P- = r _ = R_, 6_ = t?}, 

(n) 

for some constants R + , without loss of generality, and must satisfy (e.g. [10]) 

R = R+ = , [A] = 0 , [v] = 0 , [v] = 0, (12) 


where [/] is the difference of a function / at the fluid region and the value of / at the 
vacuum region evaluated on E 0 , i.e. [/] := / + |s 0 — f~ |s 0 - The first condition establishes 
that the constants R ± must agree, the condition [A] = 0 relates the constant M of the 
vacuum solution with the mass of the fluid M + (R) = M , [v\ = 0 fixes the value of the 
function at the origin and [z/] = 0 then implies that the pressure must vanish at the 
surface of the star P(R) = 0. The constant R is called a in [fy, ITU] . 

With the background configuration already matched, the first order perturbations can 
be studied. In order to cast the first order field equations in a compact form, an auxiliary 
function j is defined as 

3 ■■= e Ru+x)l2 , j~ = 1. (13) 


In the first order perturbation tensor (|2]), only the function uj(r) is involved. In fact, 
it is convenient to work with the function u(r) Q — c o(r) that satisfies the single field 
equation [5] 


1A 

r 3 dr 



+ 4#a, = o. 

dr 


(14) 
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This equation is integrated from the origin outwards, with conditions that ensure that 
the solution is regular at the origin. Thus, there is only one parameter that must be 
provided to completely determine the interior solution. In the present work it is chosen to 
be the value of Co + at the origin, denoted by u)+, although it is common in the literature 
to specify others, like the critical angular velocity for which rotational shedding occurs 


m- 


Equation (14) in the vacuum region holds for j 


1, which leads to 


to (r_) = fl 


2 J 

r _3 ’ 


(15) 


given asymptotic flatness. Choosing the spacetime gauge to first order in the fluid region 
conveniently, the functions Co + and Co~ satisfy the matching conditions [cD] = [a/] = 0 [10]. 
Hence the constant H is indeed the rotation of the fluid measured by a distant observer 
(see [TO] for a deeper discussion) . This constant along with the angular momentum of 
the star J are thus determined by 


J 


6 \dr+J R 


H — Cu + (R) + —. 


(16) 


The moment of inertia / is defined in terms of these two constants as / := J/Q. 

Regarding the second order perturbations, the change in mass involves the l = 0 

sector alone and this can be studied independently from the l = 2 sector [6j. The functions 

involved in the l = 0 sector are thus mo and ho- Although in and [10] the approaches 

taken to describe the deformation of the star are different at a fundamental level, the 

deformation is described in the three references by the value that the perturbation of the 
(2) 

pressure P ( j ; takes at the surface of the non rotating star. For convenience, this function 
is scaled as 


d 2 V+) := 


?( 2 ) I 


(17) 


2 (E(r + ) + P(r + ))' 

In the interior region the function h \is replaced by Pq 2 ' 1 and the held equations provide 


a system of first order inhomogeneous ODE’s for the set {m[J", P 0 2 }. It reads mm 


drriQ 
dr .|_ 

dP ,f 

dr .|_ 




.2 4 

r+ 


duj + 

dr+ 


2 3 ■+ 

-~3 r+J ^ 


■co 


+ 2 


Att^E + P)r + 2 ~ ( 2 ) 
r+-2M+(r+) 0 


12(r + — 2M+(r + )) 


_ m^r + 2 

(r+ -2M+(r+)) 2 1 

( + 

\ dr + J d dr + \ r. 


8nP + 



r + 3 j +2 u + ' 2 
l - 2M+(r + ) 


(18) 


(19) 


The equations must be integrated with conditions that ensure regularity at the origin. 
In addition, it is imposed that the central pressure of the non rotating configuration is 
preserved in the rotating model, so that Pq 2 ^(0) = 0 ( see [I| f° r a deeper discussion). 
The system admits the hydrostatic equilibrium first integral [6], which allows for the 
determination of . Explicitly, 

H 2) + h o ~ ^ r + 2 e~" + u +2 = 7 , 
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where 7 is a constant. Note that, once Pq 2 \o) = 0, 7 equals the value of at the origin. 
The set of functions that determine the exterior configuration to second order is 


{m 0 ,h 0 }. The held equations are (18), with j =1 and lu given by (15), and the 
following first order equation for /i ( J 


dh n 


rrin 


3 J 2 


dr_ (r_ — 2 Ad) 2 ri(r_ — 2 M) 

The asymptotically hat vacuum solution thus reads [6] 


m 0 (r_) = SM —, 


( 20 ) 


M r -) = 


SM 


+ 


J 2 


— 2M r_ 3 (r_ — 2M) ’ 


for some arbitrary constant <5M. This constant is identified as the change in mass. 

At this point the interior solution {m ( |, P 0 ' 2 ' 1 } is completely determined, while the 
exterior solution m q (20) is determined up to SM. lu order to hx it the interior and the 
exterior solutions must be related at S 0 . The full second order matching is presented in 
[ 1 10] but in here only the result regarding the functions is needed, which on Eq satisfy 


R 3 


[mo] = ~^(.R ~ 2 M)E(R)P?\R). 

This matching condition thus hxes the constant SM as 

j2 p3 

SM = m+(R) + -jp+4n—(R - 2M)E(R)P^(R). 


( 21 ) 


( 22 ) 


In order to compare the results here with [3j the change in mass is split into two contri¬ 
butions: the one given originally in reference p] and used in [3], SM ^ := m^R) + 
and the correction found in [TO] . SM ^ = 47r^(i? — 2 M)E(R)Pq 2 \R). 


3 Homogeneous stars 

When the energy density E is constant, the equations of structure (g,(@) that govern the 
(static, spherically symmetric) background configuration admit an analytical solution. In 
terms of the constant density E and the central pressure P c , that solution is given by 


P + - 
P + E 

= -A+(r+) 


P + - 

1 C ' Q 


= 1 - 


P c + E" 1 

2 M+(r+) 


8irEr+ 2 


' + 


4-7T 

M+(r+) = —Er + 3 


, v +( r +)/2 = e u +( 0)/2 


p + - 

-1 + M±mL 

Pc + E 


8-KErP 2 


-i + h±l 

Pc + E 


-1 


(23) 

(24) 

(25) 


This solution stands for the whole interior region, i.e. from r + = 0 to r + = R. The 


vacuum solution is given by (10) and extends from r + = R to infinity. The two solutions 
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are related by means of the matching conditions for the background configuration (12). 
The “continuity” of A and u' implies that M = AirER 3 /3 and P(R) = 0. 

With the background configuration already matched, it is convenient to change from 
the interior parameters {E,P C } to the exterior parameters {M,R}. Thence the solution 
(23) takes the form as given in |9, 12]. In order to present the results as in [3], the exterior 
parameters still have to be scaled with the Schwarzschild radius Rs '■= 2 M so that they 
become {Rs, R/Rs}- Inverting the relation between the parameters one finds 


E = 


8ttR 2 s 


R 

R~s 


-3 


Pr = 


8tcR 2 s 


R 

R~S 


-3 1 — 




-i 


3a 1 


1 - 


(26) 


\ R sJ 


Equation (26) implies a constraint on the background exterior parameters. In order to 
keep the central pressure finite, the following inequality must hold (3] 


0 < -M < R 0 < R < 
~ 4 


37 tE 


(27) 


In order to write the field equations for the fluid as in [3], two constants that replace 
E and R must be introduced as 


a : = 


k := 3 




8'kE' v cr 

The constant a is related to the Schwarzschild radius Rs by 


1. 


a = R 


s 


7 ? \ 3 / 2 

s 


(28) 


(29) 


Hence, given any function in units of a, it can be easily converted to units of the 
Schwarzschild radius by simply scaling it with the proper factor R/Rs- Following the 
conventions in (3], the radial coordinate in the interior region is finally substituted by 


x := 1 


i-W 


a 


(30) 


This change is well defined once the inequality (27) holds. The domain of definition of 
this coordinate is x G [0, 2/3). The radius of the star is denoted by X, which in terms of 
k is expressed as X = (2 — k)/3. Finally, in terms of x the auxiliary function j + reads 
j + = 2(1 — x)/(k + x). 

Considering the background solution (23)-(25) and the definitions introduced so far, 
the first order field equation (14) for the interior region is written in terms of the radial 
coordinate x as 

— x(2 — x){x + k) + (4x 2 — x(3 — 5k) — 5k)— -f- 4(1 + k)uj + (x) = 0. (31) 


dx 2 

The behaviour of d) + near the origin is fixed by 


u = ui’ 


1 + 4{ \ +k) x 
5k 


dx 


+ 0(x 2 


(32) 
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The second order field equations (18) and (19) are written in terms of the radial 
coordinate x as |3] 


a 


_3 duiQ 
dx 


[l - x){{2 - x)xf/ 2 f 1 _ (duP 

(.k + x ) 2 \3 [ } \dx 


8(k + 1 ) +2 
3(k + o;) 


(33) 


a 


—2 dPo 


( 2 ) 


K + 1 


da; 


(1 — o;)(k + a;) 


a -^> _ 2+(K+ „ 1)(1 3(1 tf 2 V 3 m+ (34) 


8(2 — x)x , do; + 8(k(x 1) + x) „ , 9 

a; + —--—— -co + 


(k + x)(l — a;) 2 a; 3 / 2 (2 — a;) 3 / 2 
(2 — x) 2 x 2 


3 (k + x) 2 dx 
Regularity at the origin demands [3] 


m 


a 3 u+ 2 


(x —> 0) = 


3(k + x) 3 


32y / 2(/=c + l)x 5 / 2 


3(1 — x)(k + x) 2 \ dx 


du + 


15k 3 


+ 0(x 


7/23 


P 


( 2 ) 


8a; 


a 2 ui+ 2 


lx II) = 4^ + 0(x 2 ). 


(35) 

(36) 


Once the values of the functions m) and Pq 2 ' in x — X are found, the change in 


mass is calculated using expression (22). In order to present the numerical results, it is 


convenient to express 5M divided by the mass of the background configuration M, and 
in units of J 2 /R A S . It reads 


SM J 2 


M 


= JP <2 

K s 


m r 


0 J 2 /R I) 


So 


ru , R 

+ R~S 


-3 


, R 

+ 61 s; - 1 


p, 


( 2 ) 


(,P/R%) 


(37) 


Note that the field equations for the first (31) and second order 
conditions at the origin (32) and (35), (36) are formulated for u ; + /d) 


m. 


(34) and the 
["/a 3 d)+ 2 and 

/a 2 tic 2 , and thus, depend only on one free parameter, k, from the background config¬ 
uration. This parameter k is in fact determined by the ratio of the radius of the spherical 
star to the Schwarzschild radius R/Rs by means of (28) and (29). Hence the model is 
solved just specifying a value of R/Rs■ A sequence of models with different values of 
R/Rs is explored and the results for the first and second order are presented below. 

The first order results are shown in Table [lj which includes the moment of inertia and 
the value Co + / (J/R|)|s 0 for different values of R/Rs■ For a solid sphere in the Newtonian 
regime, the normalized momentum of inertia i := I/MR 2 takes the value 2/5, which is 
achieved asymptotically (see Table [I]). The comparison with a sphere makes sense because 
the first order perturbations do not change the shape of the star. Hence, the deviation 
between this value and the values shown in Table [l] is an effect of the twisted geometry. 
These results [^] fully agree with those presented in [3j. 

The numerical results for the second order are summarized in Table [2] and Figures [l] 
and[2 In Figure^m,/ in units of J 2 /R 3 S and P^ 2 ' 1 in units of J 2 / at r + = R are shown as 
functions of R/Rs, he. /(J 2 /R 3 s )\y, q (R/R s) and P^p/ (J 2 /i?|)|s 0 (-R/ Rs) respectively. 

I 11 order not to overwhelm the notation in the subsequent discussion, let us refer to these 
two previous functions as rriQ^R/R$) and Pq^R/Rs) respectively. As mentioned in |[3J 


2 Values for R/Rs = 9/8 are not shown because although the perturbations can still be solved [3]; the 

background solution is not regular. 


































R/Rs 

/ 

i 

uj(Ro) 

1.15 

0.5105 

0.7720 

0.64391 

1.2 

0.5248 

0.7289 

0.74818 

1.3 

0.5657 

0.6695 

0.85741 

1.4 

0.6171 

0.6296 

0.89174 

1.5 

0.6758 

0.6007 

0.88714 

1.6 

0.7406 

0.5786 

0.86205 

1.7 

0.8106 

0.5610 

0.82652 

1.8 

0.8856 

0.5467 

0.78622 

1.9 

0.9653 

0.5348 

0.74441 

2.0 

1.049 

0.5247 

0.70294 

2.5 

1.534 

0.4910 

0.52376 

3.0 

2.123 

0.4717 

0.39700 

4.0 

3.604 

0.4505 

0.24623 

5.0 

5.487 

0.4390 

0.16625 

10.0 

20.91 

0.4182 

0.045821 

20.0 

81.77 

0.4088 

0.011980 

35.0 

248.1 

0.4050 

0.0039848 

50.0 

504.3 

0.4035 

0.0019668 

100.0 

2009 

0.4017 

0.00049585 


Table 1: I in units of i?|, i := I /MR 2 and lj(Rq) / (J/ R 3 S ) for some values of R/Rs- 


these are not monotonic functions and they both present a maximum, the first one at 
R/Rs rv_/ 1.29 and the second at R/Rs ~ 1.82. Note that the function is negative for 
small values of R/Rs and this implies that the average deformation of the star to second 
order |^] can be either negative or positive, so that the star may show either contraction or 
expansion to second order depending on the background parameters. It is worth noting 
that m ot, and P^ attain values of the same order and the ratio is about 1.63 for big values 
of R/Rs ■ For a more detailed discussion we refer the reader to the original work PJ. 

Finally, the results including the corrected change in mass and their comparison with 
those presented in [3] are shown in Table [2] and Figure [2j In Table [2] some values of 
the change in mass as a function of R/Rs are presented. In the second column the 
value of SM^/M, which corresponds to the 5M/M given in [3], is shown, whereas the 


third column includes the correct change in mass (22) and, lastly, the fourth column 


shows the fraction of the change in mass that corresponds to the correction. In fact, the 
correction becomes the dominant contribution in 5M as the quotient R/Rs increases. The 
behaviour of 5M/M is shown for a wide range of the variable R/Rs in Figure [2] SM/M 
presents a maximum at R/Rs ~ 2.81 and then decreases more slowly than the SM^/M 
presented in [3], which decays monotonically. The amended 8M/M and SM^/M agree 
for R/R s ~ 1.27158. 

There is a combination of two facts that makes the correction 8A / R C ' 1 not only notice¬ 
able, but also dominant in homogeneous stars. On the one hand, as shown in Figure [lj 
/riot, and P 0 t| are quantities of the same order. On the other hand, in formula (37), the 


coefficient with -Foil scales linearly with R/Rs- These results definitely reveal that the 


? The relation of the pressure and the shape of the star is addressed in |)j. 
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R/Rs 

SM^/M 

5M/M 

\8MW\/8M 

1.15 

3.454 

2.348 

0.4711 

1.2 

3.412 

2.725 

0.2524 

1.3 

3.225 

3.506 

0.0803 

1.4 

2.993 

4.246 

0.2952 

1.5 

2.757 

4.904 

0.4379 

1.6 

2.533 

5.474 

0.5372 

1.7 

2.327 

5.954 

0.6091 

1.8 

2.140 

6.355 

0.6632 

1.9 

1.971 

6.684 

0.7051 

2.0 

1.819 

6.951 

0.7383 

2.5 

1.259 

7.631 

0.8350 

3.0 

0.9163 

7.690 

0.8808 

4.0 

0.5440 

7.173 

0.9242 

5.0 

0.3588 

6.482 

0.9446 

10.0 

0.09493 

4.065 

0.9766 

20.0 

0.02437 

2.259 

0.9892 

35.0 

0.008049 

1.349 

0.9940 

50.0 

0.003966 

0.9608 

0.9959 

100.0 

0.001008 

0.4901 

0.9979 


Table 2: The change in mass SM^/M typical from the literature, the amended change 
in mass 5M/M and the fraction of the correction with respect to the total change in mass 
are presented for different values of R/Rs- 


correction to the change in mass is important in homogeneous stars and it can not be 
neglected at all. Again, the correction to 5M depends crucially on the value of the energy 
density at the surface of the static and spherically symmetric configuration, and hence, it 
does not contribute, in principle, to SM in the usual EOS’s for neutron stars. However 
actual models of neutron stars show a residual value associated with the numerical error 
present when the radius of the star is determined by P(R) = 0. Therefore, the correction 
to the change in mass might be relevant in order to improve the precision in the determi¬ 
nation of the mass of the rotating configuration. This fact regarding the usual EOS’s for 
neutron stars and other cases of interest are currently under study. 
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Figure 1: The perturbation to the pressure to second order Pq^R/Rs) and m^(R/ Rs)- 



Figure 2: The original and the amended changes in mass versus the normalized radius of 
the static star R/Rs- 
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